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Preface 


"This  is  Blue  Squadron  Leader.  They're  jamming 
our  communications  over  the  target  — switch 
to  SS  mode." 

This  could  have  been  a line  from  Star  Wars  or  a mili- 
tary pilot  in  our  next  armed  conflict.  Times  have  long  since 
past  when  wars  were  won  by  the  side  with  the  bigger  stick. 

In  this  age  of  electronic  wizardry,  the  side  that  pushes  the 
state-of-the-art  has  the  upper  hand.  This  paper  investigates 
the  use  of  charge-coupled  devices  (CCDs)  in  spread  spectrum 
receivers  which  even  up  the  size  of  the  sticks  in  the  jamming 
environment.  Since  all  communications  texts  use  different 
notation,  some  familiar  benchmark  results  are  presented  to 
provide  a point  of  calibration. 

Working  on  this  thesis  with  the  assistance  of  Captain 
Stan  Robinson  as  advisor  was  a refreshing  experience.  It 
would  not  be  fair  to  thank  any  one  person  for  the  behind  the 
scenes  work  in  putting  this  thesis  together;  so  thanks  are 
due  to  Mrs.  Lana  Apana  for  making  these  words  and  equations 
come  to  life.  I am  endeared  to  Susan,  Michael  and  Cory  for 
their  total  support  during  my  entire  AFIT  program. 
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Abstract 


A phase  comparison  receiver  is  implemented  at  baseband 
with  charge-coupled  device  (CCD)  correlators  as  matched  fil- 
ters. The  received  signal  is  a pseudo-noise  (PN)  phase  coded 
spread  spectrum  waveform.  The  effects  of  doppler-induced 
frequency  instabilities  on  the  probability  of  error  perform- 
ance is  examined  with  the  differential  phase  modeled  as  a 
Gaussian  process.  Numerical  results  are  presented  and  when 
the  variance  is  greater  than  0.1  radians2  the  performance 
degrades  rapidly.  An  analytic  approximation  is  given  for  prob- 
ability of  error  for  specific  cases  and  is  useful  for  small 
values  of  variance  and  low  signal-to-noise  ratios.  In  the 
limit  when  the  variance  approaches  zero,  the  receiver  perform- 
ance expression  is  equivalent  to  the  differentially  coherent 
phase  shift  keyed  (DPSK)  receiver.  If  the  reference  is  perfect 
in  addition  to  the  zero  phase  variance  the  performance  is  the 
same  as  a phase  coherent  receiver  with  a perfect  reference. 


I.  Introduction 

< 

Electronic  warfare  has  come  a long  way  since  the  days 
when  the  bad  guys  would  cut  the  telegraph  lines  to  the  unfor- 
tunate town.  Electronic  warfare  cannot  be  ignored;  it  is  a 
significant  present-day  reality  which  must  be  countered  with 
state-of-the-art  technology.  The  simplest  method  of  denying 
our  pilots  the  use  of  their  radio  communications  is  to  over- 
power (jam)  their  transmissions.  One  method  of  increasing 
the  survivability  of  radio  communications  in  the  jamming 
environment  is  the  spread  spectrum  modulation  technique. 

Spread  spectrum  (SS)  systems  were  built  as  early  as  the  late 
1940s  (Dixon,  1976:6-7).  The  hardware  complexity  of  these 
early  receivers  filled  entire  rooms  with  radio  equipment  limit- 
^ ing  the  practicality  of  their  use  in  many  applications.  Device 
technology  has  advanced  since  then  to  the  point  where  spread 
spectrum  transceivers  can  be  contained  in  handheld  units  (Dixon, 
1976:10)  . 

One  particular  new  device  that  will  be  useful  in  future 
receiver  design  is  the  charge-coupled  device  (CCD) . The  CCD 
was  invented  in  1970  (Melen,  1977:1),  but  the  need  for  such  a 
device  was  overdue  and  commercially  available  models  are  now  in 
production.  Unlike  a number  of  spread  spectrum  correlation 
techniques,  the  CCD  operates  in  the  baseband  frequency  range  of 
1-10  MHz  rather  than  at  RF  or  IF  frequencies.  Baseband  process- 
ing of  coherent  signals  requires  that  either  the  RF  carrier  be 
phase- locked  to  the  local  oscillator  or  that  both  quadratures  of 
( of  the  RF  signal  be  processed  and  recombined  at  baseband  to 
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eliminate  the  phase  differences  between  the  local  oscillator 
and  received  signal.  In  this  paper,  the  phase  compensation 
will  be  accomplished  using  the  latter  method  where  both  quad- 
ratures will  be  used  for  baseband  processing. 

The  phase  offset  between  the  received  signal  and  receiver 
local  oscillator  is  assumed  to  be  a result  of  air-to-air  or 
air-to-ground  communication  and  is,  therefore,  referred  to  as 
doppler  shift.  However,  the  model  for  phase  instability  can 
be  extended  to  analyze  the  effects  of  phase  shift  due  to 
oscillator  drift,  atmosphere  and  other  common  disturbances. 

In  this  paper,  spread  spectrum  receivers  using  CCD  code  correla- 
tors are  developed  and  analyzed  for  binary  signalling  probability 
of  error  performance. 

Background 

There  are  four  basic  techniques  for  spread  spectrum  modula- 
tion discussed  in  the  literature  (Dixon,  1976) : direct  modulated 
sequences,  frequency  hopping,  pulsed-FM  or  chirp,  and  time  hopping. 
The  receivers  discussed  in  this  paper  use  exclusively  the  direct 
sequence  format.  One  method  of  generating  a direct  sequence  spec- 
trum from  binary  data  is  to  multiply  a data  bit  by  a pseudo-random 
code.  The  code  bits  or  chips,  as  they  are  more  commonly  known, 
are  of  much  shorter  duration  than  the  data  bits.  The  resultant 
signal  is  used  to  modulate  an  RF  carrier.  Since  the  chip  time 
dominates  or  masks  the  data,  the  side  lobes  of  the  frequency 
spectrum  are  spaced  farther  apart.  Additionally,  the  power  in  the 
signal  remains  the  same  as  the  unspread  signal;  hence,  the  power 
( density  is  reduced.  Figure  1 shows  the  frequency  spectrum  of  a 
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Figure  1.  Data  Modulated  Carrier  and  PN  Modulated  Carrier 

direct  sequence  modulated  carrier  versus  the  data  modulated 
carrier. 

The  price  for  increased  resistance  to  jamming  is  an 
increase  in  the  complexity  of  processing  requirements.  The 
CCD  has  many  desirable  features  necessary  in  the  processing 
of  spread  spectrum  signals.  The  CCD  provides  a flexible  inter- 
face between  continuous  analog  signals  and  sampled  data  process 
ing.  Additionally,  it  is  easily  adapted  to  microprocessor  con- 
trol which  adds  flexibility  in  adaptive  processing  and  decoding. 

o 
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Objectives 


The  scope  of  this  paper  includes  design  and  (binary) 
probability  of  error  performance  of  a spread  spectrum 
receiver  using  CCD  code  correlators.  The  binary  probability 
error  performance  is  considered  for  different  energy  per 
bit-to-noise  spectral  density  ratios  and  doppler  rates.  The 
receiver  noise  is  white  Gaussian  noise  (WGN)  and  the  doppler 
is  due  to  the  relative  motion  between  receiver  and  trans- 
mitter. 

The  receivers  analyzed  do  all  processing,  including 
doppler  phase  tracking,  at  baseband.  The  local  oscillator/ 
mixer  will  always  be  tuned  to  the  transmit  frequency.  How- 
ever, it  is  not  a voltage  controlled  oscillator  nor  is  it 
phase-locked  via  a phase-lock  loop  (PLL) . 

Approach 

This  report  is  organized  into  three  main  sections  with 
a chapter  to  cover  each  topic.  Chapter  II  explains  the  CCD 
correlator  and  the  baseband  output  signal  that  is  common  to 
all  receivers.  A set  of  benchmark  receivers  are  discussed  in 
Chapter  III  as  a basis  of  comparison  for  error  performance. 

The  spread  spectrum  receivers  are  also  developed  and  analyzed 
for  error  performance  in  Chapter  III.  Conclusions  and  recom- 
mendations are  included  in  Chapter  IV.  The  appendices  contain 
some  of  the  more  involved  and  detailed  derivations  necessary  for 
analysis.  When  the  results  of  a calculation  or  derivation  are 
included  in  the  appendix,  it  is  mentioned  in  the  corresponding 
section  of  the  text. 
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Chapter  II.  Received  Signal  and  CCD  Model 


The  received  signal  and  the  CCD  code  correlators  are  dis- 
cussed in  this  chapter.  Each  signal  quadrature  is  processed 
through  the  code  correlators,  with  the  resultant  sampled  data 
quadrature  values  processed  by  the  baseband  receiver. 

Received  Signal  Model 

The  received  waveform  consists  of  a binary  phase  shift 
keyed  (BPSK)  signal  component  and  a white  Gaussian  noise  (WGN) 
component  due  to  receiver  noise  and  broadband  jammers.  The 
received  waveform  is  given  by  Eq  (1) . 

r(t)  = s(t)  + n(t)  oitlTb  (1) 


where 

s (t)  is  the  signal 

n(t)  is  the  noise 

s 

The  noise  is  WGN  with  a double  sided  spectral  density  of 
1/2  Nq  watts/Hz.  In  addition  to  natural  sources  of  noise, 
additional  noise  power  from  a jammer,  assuming  the  jammer  is 
using  broadband  noise  is  included  in  the  1/2  No  term.  The  sig- 
nal portion  of  the  received  signal  is  given  by: 

s(t)  = /2Ad(t)c(t)  cos  (2irf0 1 + 0 (t)  ) (2) 

where 
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A is  the  amplitude  of  the  received  signal 
d(t)  is  the  data  waveform,  di_  (t)  = -d2  (t)  - 1 
c(t)  is  the  PN  chip  waveform,  |c(t)|  = 1 

0 is  the  random  phase  difference  between  the  received 
signal  and  local  oscillator  due  to  transmit  and  receive  oscil-  . 
lator  phase  difference,  doppler,  and  delay 

Each  term  is  assumed  to  have  the  following  properties 

throughout  the  analysis.  The  amplitude  term  A is  a constant. 

The  data  waveform  d(t)  is  equiprobable  and  has  a value  di (t)  = 

-d2 (t)  and  |d(t)|  = 1 for  a bit  time  Tfa  seconds.  Similarly, 

the  chip  waveform  is  c(t)  with  |c(t)|  = 1 for  Tfa  seconds 

and  varies  at  a chip  rate  R which  is  much  faster  than  the  bit 

c 

rate  R^  = 1/T^  . The  carrier  frequency  fo  is  constant  and 
known  to  the  receiver.  The  description  of  the  phase  0 ( t ) is 
covered  in  the  next  section. 

Statistical  Model  for  Phase  Offset 

The  change  in  phase  0 (t)  is  assumed  to  be  due  to  the 
relative  motion  between  the  transmitter  and  receiver.  For  example, 
if  the  receiver  and  transmitter  are  aboard  two  different  aircraft, 
the  relative  velocity  between  the  two  aircraft  depends  on  their 
speed  and  direction.  If  the  aircraft  are  traveling  alongside 
each  other,  their  relative  velocity  may  be  zero,  in  which  case 
the  instantaneous  doppler  frequency  is  d0(t)/dt  =0  Hz.  If 
the  aircraft  are  traveling  away  from  or  toward  each  other  along 
the  same  velocity  vector,  the  instantaneous  frequency  is  given  by: 
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de (t)  = 


fo  Hz 


(3) 


where 

Vi  and  V2  are  the  velocity  vectors  of  aircraft  1 and 
2 respectively 

c is  the  speed  of  light 
f is  the  carrier  frequency 


A statistical  model  for  the  phase  process  is  assumed  mod- 
erately wide  sense  stationary  over  several  bit  times  such  that: 


( 


E [0  ( t ) ] = 0o  radians  (4) 

where 

E[xJ  is  the  expected  value  or  ensemble  average  of  the 
waveform  x(t)  (Papoulis,  1965:138-151) . 


The  generic  receiver  is  a phase  comparison  receiver  where 
the  previously  received  phase  is  used  as  the  reference  of  the 
present  phase,  A differential  phase  process  generated  by  sampling 
0 (t)  at  the  bit  times  is  of  primary  interest  rather  than  the 
density  of  0 (t)  . The  differential  phase  process  A0 (t)  will 

be  defined  in  terms  of  the  known  parameters  of  f^  and 


i 


4 


( 


7 


the  correlation  properties  of  0 (t)  . The  fluctuation  in  6 (t) 

is  assumed  to  be  due  only  to  doppler  as  a result  of  relative 
vehicle  velocity.  A new  zero  mean  process  e ^ C t ) thus  is  given 
by: 


ed(t)  = e (t)  - e0 


(5) 


The  variance  of  the  first  derivative  will  be  defined  as 
the  maximum  mean  squared  doppler  frequency  (Papoulis,  1965: 
339)  : 


(6) 

(7) 


where 


-82Red(T) 

3T2 


T = 0 


T = |t-t| 


max*  f o 


(8) 


x 


Rx(T)  is  the  correlation  of  the  function  of  the  process 
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The  correlation  function  of  0.(t)  is  assumed  to  be  a 

a 

Gaussian  shaped  function  as  the  second  derivative  exists  every- 
where and  is  given  by: 


R (T)  = a 2 exp  (-T2/n2) 
0d  6d 


(9) 


The  value  of  a 2 can  be  defined  in  terms  of  the  second 
0d 

partial  derivative  of  R0<^(T).  On  computing  the  second  partial 
as  in  Eg  (7)  and  solving  for  a 2 the  result  is: 


°ed2  " f <^2 


(10) 


The  correlation  time  is  given  by  (Lindsey,  1973:313) 


Tc  °0. 


2 - 


I Rea<T>  | 


(ID 


Using  Eq  (9)  and  solving  for  T gives: 

c 


T = — /iT 

c 2 


(12) 


The  CCD  code  correlators  require  the  phase  to  be  relative- 
ly constant  over  a symbol  interval  Tfa  a point  that  will  become 
apparent  in  the  next  section  on  the  CCD  code  correlators.  Typi- 
cally, the  phase  can  be  considered  constant  over  a symbol  inter- 
val if  (Lindsey,  1973:314): 


T 


c 


> 


4 Tt 


(13) 
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In  terms  of  n 


0 < T.  < - /H  = T c (14) 

_ b-  8 — 

The  value  of  is  limited  by  Tc  which  is  in  turn 

a function  of  the  dynamics  of  the  aircraft  in  relation  to  each 
other  during  a typical  mission.  Useful  values  of  Tc  are 
not  presently  available  and  are  left  unspecified  as  is  f/d. 

The  results  are  then  applicable  in  the  most  general  case.  The 
differential  process  is  given  by: 

A6  (t)  = 0 ( t)  - 0(t  - Tb)  (15) 

The  statistics  of  A0  are 

E [ A 0 ] = 0 (16) 

and  (Papoulis,  1965:336-381): 

“U1  ■ 2Rea(T)  - Red(lTb‘Tl)  ' Rea<  Tb+T  > <17) 

Since  A0  is  a zero  mean,  the  variance  is  given  by: 

°A02  = RA0(o)  = 2o0d  (1  - exp(-Tb2/n2))  (18) 

Using  Eq  (14) , the  admissible  range  of  the  variance  of 
A0  is  thus  given  by: 
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0 < a 


- A0  - 


< .la 


0d 


2 2 

2 = f . T_ 


(19) 


5ir 


The  range  of  a.  2 is  now  defined  in  terms  of  measur- 
able  parameters  Tc  and  f^.  Information  on  a useful  density 
of  A0  was  not  available  at  this  writing;  however,  it  is 
reasonable  to  assume  that  the  A0  is  distributed  as  a zero 
mean  Gaussian  random  variable  with  aAg2  given  by  Eq  (18). 

A useful  density  which  exhibits  properties  similar  to  a Gaus- 
sian  density  but  is  limited  to  + tt  is  given  by: 


*««»  - exs,1;0cu » - 


(20) 


where 


I is  modified  Bessel  function  of  the  first  kind  of 
0 


order  zero 


a is  approximately 


for  a > 4 or  a192  < _25 


A0 


(Viterbi,  1966:92-95) 

The  density  in  Eq  (20)  will  be  used  for  the  statistical 
model  of  A0  as  in  addition  to  its  Gauss-like  characteristics, 
it  has  well  defined  higher  order  moment  properties  in  cos0  and 
sin0  which  will  be  of  value  when  investigating  the  quadrature 
representations  of  the  received  signal  (Lindsey,  1973:36). 

CCD  Model 

CCDs  have  only  recently  become  available  to  designers 
in  commercial  quantities  in  both  analog  and  digital  configurations 
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(Melen,  1977) . A wide  variety  of  CCD  models  exist  for  analog 
signal  processing,  consequently  the  mathematical  model  derived 
for  analysis  is  based  on  a specific  correlator;  the  General 
Electric  CR  series  of  CCD  correlators  (Buss,  1973) . 

A CCD  configured  as  an  analog  correlator  consists  of  a 
number  of  charge  storage  bins  or  stages.  The  amount  of  charge 
in  each  bin  is  proportional  to  the  input  signal  level  and  an 
integration  time  during  which  the  bin  is  filled  with  a charge. 

A typical  CR  correlator  consists  of  32  stages  and  several  CRs 
can  be  cascaded  to  form  correlators  of  multiples  of  32  stages. 

The  CR  correlators  operate  on  a moving  reference  rather 
than  a moving  charge  principle.  The  moving  reference  CCD  con- 
figuration reduces  the  number  of  singal  packet  shifts  which 
improves  efficiency  of  the  device  (Buss,  1973:83-90).  A block 
diagram  of  the  correlator  is  shown  in  Figure  2. 

In  order  to  simplify  the  analysis,  th.e  effects  of  atten- 
tuation  and  noise  added  by  the  CCD  will  not  be  included  in  the 
derivation  of  the  output  statistics  of  the  code  correlators. 

The  inclusion  of  the  effects  of  attentuation  and  noise  in  the 
CCD  would  not  drastically  change  any  results  as  the  loss  through 
a CCD  of  32  stages  is  on  the  order  of  fractions  of  a dB  and  the 
CCD  internal  noise  is  considerably  less  than  the  receiver  input 
noise  (Buss,  1973:90). 

The  correlators  are  connected  to  form  the  desired  receiver 
as  shown  in  Figure  3.  The  output  of  the  multipliers  with  s(t) 
as  the  input  ignoring  double  frequency  terms  (as  the  frequency 
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Figure  2.  Moving  Reference  CCD  Code  Correlator. 

characteristics  of  most  multipliers  do  not  cover  such  a wide 
frequency  range)  is: 

s(t)  (/I  cos  2ir f t)  = S_(t)  (21) 

o J- 


then: 


Sjtt)  = Ad(t)c(t)  cos  0(t) 


similarly 


Sg(t)  =*  Ad(t)c(t)  sin  0(t) 


The  code  chips  c(i)  are  + 1 and  assumed  to  be  properly 


1 


similarly 


Q (k)  = ANT  d(k)  sin  6 (k) 
s c 


(26) 


where 


N is  the  number  of  correlation  stages  (also  the  num- 
ber of  chips  in  code  bit) 

(k)  is  the  notation  for  (kT^)  and  the  is  omitted  to 
simplify  notation 


Because  of  the  restrictions  imposed  on  the  phase  process 
8(t),  the  phase  is  constant  over  a symbol  interval  of  seconds 

and  cos  0 (k)  is  given  by  (Lindsey,  1973:312): 


cos  0 (k) 


k + 1)  T. 
cos  0 (t 


kTb 


)dt 


(27) 


with  a similar  expression  for  sin  0 (k) . 

The  statistic,  .f  the  output  noise  Nj (k)  are: 


E [Nt (k) ] = E 
1 \ 


r n r 

X L V^/n 


'•T.  + T 
b c 


(t)  cos  2irf Q + dt 


0 k Tb  + (i-l)Tc 


(28) 


= 0 


where 
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(32) 


SNR 

c 


E[I(k)  2] 
EINj.  (k) 2] 


SNR 

S 


E [Q  (k) 2 ] 
E[NQ(k)  2] 


(33) 


The  signal  to  noise  input  to  the  received  signal  processor  is 
defined  as  the  SNRC  with  0(t)  =0.  The  output  of  each  quad- 
rature is  given  by: 


I (k)  = Is(k)  + NI(k) 


(34) 


Q (k)  = Qs(k)  + NQ(k) 


(35) 


Note  that  the  expected  value  of  I(k)  and  Q(k)  in  Eqs  (32)  and 
(33)  is  the  signal  portion  of  Eqs  (34)  and  (35),  respectively. 

Using  the  value  of  Ig(k)  from  Eq  (25)  with  0 (k)  = 0 
and  noting  that  d(k)2  = 1 the  signal  to  noise  is: 


SNR 


A2N2T  2 
c 

N N T 
o c 


(36) 


Recalling  that  NTC  = Tfa  , we  will  see  later  that  the 
performance  will  be  expressed  in  terms  of  energy -per -bit  to 
noise  spectral  density  given  by: 


R 


A2NTt 


N. 


(37) 
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The  model  for  phase  and  the  expression  for  energy  per 
bit-to-noise  is  used  throughout  the  remainder  of  this  paper. 
The  output  of  the  correlators  in  Eqs  (34)  and  (35)  is  the 


input  signal  used  by  the  receivers  throughout  the  remainder 
of  this  paper. 


Chapter  III.  Receiver  Design  and  Performance 


The  spread  spectrum  receivers  are  designed  and  analyzed 
in  this  chapter.  Prior  to  embarking  on  the  development  and 
performance  of  specific  receivers,  the  results  for  some  stan- 
dard benchmark  receivers  are  presented  to  provide  a point  of 
reference.  The  approach  to  the  design  of  the  spread  spectrum 
receivers  will  begin  with  the  development  of  a generic  receiver 
with  generalized  performance  results.  The  generic  model  is 
followed  by  some  special  cases. 

Benchmark  Receivers 

The  results  for  two  benchmark  receivers  are  presented 
in  this  section.  The  first  is  a coherent  receiver  with  a per- 
fect reference,  and  the  second  is  a differentially  phase  shift 
keyed  (DPSK)  receiver  with  the  signal  phase  constant  during  a 
symbol  interval  of  2T^  seconds. 

Coherent  Benchmark  Receiver 

The  classical  phase  coherent  receiver  consists  of  a local 
oscillator  tuned  to  the  received  signal  carrier  with  no  phase 
offset.  A block  diagram  of  the  receiver  is  shown  in  Figure  4, 
where  r(t)  is  given  by  Eq  (1).  As  before,  the  clock  is  syn- 
chronized to  integrate  each  chip  for  a period  of  Tc  seconds 
and  the  data  bit  is  correlated  and  sampled  every  T^  seconds. 
The  output  I(k)  is  given  by: 


I(k) 


Ad  (k)  NT  + NT(k) 
c I 


(38) 
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Figure  4.  Coherent  Benchmark  Receiver. 


Using  the  results  from  Chapter  II  and  from  Figure  4 
6 (k+i)  = 0 the  performance  is  given  by  (Van  Trees,  1968: 
98-102) : 


Pe  = 1/2  erfc  (/r~ ) 


where 


R is  defined  in  Eq  (37)  as  only  the  SNR  tern  is 


present  and 


erfc (y) 


co 

■ / 


exp  [x2]  dx 


The  next  logical  step  is  to  graph  Pg  versus  R, 
but  this  step  will  be  postponed  until  both  benchmark  results 
are  presented. 

DPSK  Benchmark  Receiver 
1 

The  DPSK  receiver  uses  the  previous  bit  as  a reference 
for  the  next  bit.  The  assumption  is  made  that  for  this 
receiver  s(t)  is  differentially  encoded  at  the  transmitter. 
Since  the  previous  data  bit  is  used  for  a reference,  a locally 
generated  synchronized  carrier  is  not  necessary,  thus  coherent 
carrier  detection  synchronization  circuits  such  as  phase-locked 
loops  are  not  necessary.  The  diagram  for  the  baseband  DPSK 
receiver  is  shown  in  Figure  5. 

The  performance  for  the  DPSK  receiver  is  derived  in  many 
texts  (Downing,  1964:184-187);  (Lindsey,  1973:245-252);  and 
(Spilker,  1977:336).  In  all  of  these  texts,  the  assumption  is 
made  that  the  phase  is  constant,  but  unknown  over  2Tj3  seconds 
and  the  only  corruption  of  phase  is  due  to  Gaussian  noise.  In 
terms  of  the  model  for  the  doppler  phase  given  in  Chapter  II, 
the  variance  °A02  = 0.  Various  methods  are  given  for  deriving 
the  performance  and  the  result  in  all  cases  is: 

Pe  = 1/2  exp  (-R)  (40) 

The  Pe  versus  R is  plotted  for  both  benchmark  receiv- 
ers in  Figure  6. 

( 
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Figure  5.  DPSK  Benchmark  Receiver. 

Generic  Spread  Spectrum  Receiver 

A generalized  diagram  of  the  spread  spectrum  receivers 
analyzed  in  this  section  is  given  in  Figure  7.  An  expression 
for  performance  will  be  presented  in  general  and  special  cases 
with  typical  values  graphed  to  compare  performance. 

The  decision  variable  W(k)  in  Figure  7 can  be  defined 
as  consisting  of  two  possible  signals.  The  signals  are  a 
reference  vector  composed  of  $(k)  and  Q(k)  and  a data  meas- 
urement vector  l(k)  and  Q(k)  in  vector  notation  rewritten 
as: 
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The  optimum  decision  is  given  by  (Wozencraft,  1968:214- 


217)  : 


T 

U(k) 


_ 1 
S0 (k)  > 

0 


(43) 


where 


x ^ y is  the  double  inequality  when  x>y  a is  chosen 
< 
b 

and  when  x<y  b is  chosen 

x Ty  is  the  inner  produce  of  vectors  x and  y 

The  characteristics  of  U depend  on  the  configuration 
of  the  particular  receiver.  In  each  of  the  special  cases  of 
receivers  developed,  the  form  of  the  vector  U will  change 
and  So  will  be  defined  in  Eq  (42)  . The  characteristics  of 

A A 

I and  Q will  be  altered  to  develop  various  forms  of  the  vector 
U . 

The  error  performance  of  the  generic  receiver  is  given 
by  (Stein,  1964:43-51): 


peUe  = Q(/i'  /E”  °-ll- 


ozs  + ozu 


exp 


a+b 


I0 (/ab)  (44) 


( 


with  an  alternate  form  given  by: 


PelA6  = 1/2  ^ + Ql/*' 


■?2s  " °2° exp  [-  2+k  ] i (/ib) 

2 (oj  + o2u)  2 


(45) 


where 


( 


Q is  the  Marcum  Q function  and  is  discussed  in  Appen- 
dix A. 

I0  is  the  modified  Bessel  function  of  the  first  kind  of 
order  zero  and  is  also  discussed  in  Appendix  A. 

o2s  is  the  variance  of  the  data  vector  terms  I (k)  and 


Q(k) 


r2u  is  the  variance  of  the  reference  vector  terms  ^(k) 


and  Q(k) 


The  terms  a and  b are  given  in  general  as  (Stein, 
1964:50) : 


E [I  (k)  ] 2 -i-  E [Q  (k)  ] 2 + E [?  (k)  ] 2 + E [Q(k)  ] 2 


u 


2 

°s  °u  L 


E (I  (k)  ] E[I(k)]  + E (Q  (k)  ] E r6  (k)  ] 


'] 


(46) 


E [I  (k)  ] 2 + E [Q  (k)  ] 2 + E [3  (k)  ] 2 + E[ft(k)]2  * 


u 


°s°u 


[e  [I  (k)  ] E [I  (k)  ] + E [Q  (k)  ] E [$  (k)  ] ] | (47) 
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where 


E [ • ] is  the  average  over  the  Gaussian  noise  statistics 
conditioned  on  the  doppler  phase  (A0)  known 

r 

The  terms  a and  b as  defined  in  Eqs  (46)  and  (47) 
assume  the  doppler  phase  A0  in  Eqs  , (25)  and  (26)  to  be 
some  constant  value;  hence,  the  phase  is  not  averaged  in  the 
expected  value  in  Eqs  (46)  and  (47) . The  terms  a and  b 
can  be  written  in  terms  of  R as  given  in  Eq  (37)  for  SQ 
and  0R  for  the  energy  per  bit-to-noise  spectral  density  of 
U.  The  terms  a and  b with  A0  known  are  rewritten: 


0) 

II 

MW 

[ (1  + 6)  - 2/ff  cos  A 0 ] 

(48) 

K - R 
b ” 2 

(1  + 6)  + 2/ff  cos  A0  ] 

(49) 

where 


A0  is  the  angle  between  the  signal  components  of  Sq 

and  U 

6 is  the  ratio  of  energy  per  bit-to-noise  spectral  den- 
sity ratio  of  the  U reference  and  S data  vector. 

o 

The  term  6 is  considered  to  be  greater  than  or  equal 
to  one  for  all  cases  considered.  It  is  tempting  at  this  point 
in  the  analysis  to  average  cosA0  over  the  statistics  for 
A9  in  Eq  (20);  however,  the  Q function  is  not  a linear 
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operation,  and  the  expected  value  operation  on  A0  cannot 
be  interchanged.  The  required  operation  for  average  probabil- 
ity of  error  performance  is  given  by: 


pe  PeU6  £<a9,d4e 


(50) 


where 


PelA6  con^^ti°na^  probability  of  error  expression 

in  Eq  (44)  or  (45)  with  A0  known 

The  density  f(A0)  is  rewritten  from  Chapter  II  Eq  (20) 
for  reference: 


< 


( 


f (A0)  = g?4?cosAA)  (51) 

/ttIq  la) 

The  average  given  on  the  right  side  of  Eq  (50)  is  not 

a trivial  operation  and  a closed  form  solution  for  the  general 

case  possibly  does  not  exist.  There  are  a few  limiting  cases 

and  special  configurations  of  the  generic  receiver  where 

approximations  for  error  performance  are  useful.  A realistic 

operational  range  of  interest  is  the  required  energy  per  bit- 

to-noise  spectral  density  that  has  an  average  probability  of 

_3  _ 5 

error  on  the  order  of  10  or  10  bits  per  error. 

The  terms  a and  b were  written  with  the  parameter 
B as  a ratio  between  the  signal  bit  and  reference  bit.  It 
is  conceivable  that  the  filter  in  Figure  7 could  smooth  the 
previously  received  vectors  and  reduce  the  variance  of  the 
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(k-1) 


(k-1) 


U vector  in  which  case  3 would  be  greater  than  one.  In 

the  case  where  the  filters  are  absent  and  there  is  no  data 

feedback,  the  typical  DPSK  receiver  results.  The  term  a 

in  Eq  (51)  is  also  a variable  parameter  in  system  performance 

controlled  by  the  coherence  time  of  the  phase  A0(t)  and  the 

data  time  T^  . If  the  filters  in  Figure  7 could  predict 

the  phase  of  SQ  , then  it  is  conceivable  that  the  variance 

on  the  phase  could  be  reduced,  thus  increasing  a ^ since  a = 

— \ . The  larger  the  term  a is,  the  greater  the  probabil- 
a 2 / 

ity  that  A 6 is  close  to  zero  radians.  The  a and  6 
terms  can  conceptually  be  varied  independently  to  determine 
which  parameter  has  the  greatest  effect  in  improving  error 
performance  without  specifying  the  filters  in  Figure  7. 

In  the  next  section,  various  configurations  of  the  gen- 
eric receiver  of  Figure  7 will  be  considered  and  analyzed  for 
error  performance.  In  the  receivers  where  data  feedback  is 
present,  the  receiver  is  assumed  to  be  operating  with  perfect 
data  feedback.  Additionally,  the  variance  of  the  vector 
components  SQ  and  U will  be  assumed  equal  although  0 
may  be  greater  than  one  (Stein,  1965:45);  with  this  assump- 
tion, only  the  first  half  of  Eq  (45)  is  necessary  for  error 
performance,  as  o2  - oo2  = 0. 

Generalized  Performance  Results.  A closed  form  solu- 
tion to  Eq  (50)  possibly  does  not  exist  in  which  case  the 
only  other  alternatives  are  bounds,  approximations,  or  numer- 
ical integration.  The  latter  alternative  will  be  presented 
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first  to  provide  a point  of  comparison  for  an  approximation 
which  is  useful  in  a limited  region. 

The  Q function  routine  was  not  available  in  the  com- 
puter library  functions;  however,  a simple  recursive  method 
used  to  generate  the  Q function  tables  commonly  used  in 
radar-detection  theory  (Brennan,  1965:312-313)  was  easily 
developed  using  a FORTRAN  program  and  the  recursion  relation 
given  in  Appendix  B. 

The  probability  of  error  was  calculated  using: 

Pe  = 1/2  [ 1 -/**  Q(/b,  /a)  f (A0)  dA6 
J-  7T 

+J  * Q(/a,  /b)f(A0)dA0  j (52) 

The  numerical  integration  technique  requires  that  the 
integrand  be  a function  of  one  variable  and  the  functions  in 
the  integral  are  rewritten: 


F (x)  = Q(  /R/ZC(l+0)  + 2^-  cos 

/572-TTTfeV  TJn  cos  4 1?  X>  <»l 

The  probability  of  error  in  Eq  (52)  can  be  rewritten 
in  terms  of  one  variable  of  integration  as 


Pe  = 1/2 


[l  -f'  F+  (x)  dx 

y[  * F"  (x)  dx  1 

[ J 

1 

(54) 
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where 


F+(x)  is  given  in  Eq  (53)  with  the  plus  in  the  first 
radical  and  a minus  in  the  second  radical  in  Eq  (53) . 

F (x)  is  the  opposite  sign  under  the  radical  in 
Eq  (53)  as  F+(x) 


The  only  variable  of  integration  in  Eq  (54)  is  x and 
F (x)  can  be  written  as  a function  subroutine  suitable  for 
use  in  a numerical  integration  library  subroutine.  Unfor- 
tunately, with  even  a large  computer  such  as  the  CDC  6600, 

each  integration  takes  approximately  .3  seconds  (with  = 

_1  5 

10  ; see  Appendix  B)  which  takes  60  seconds  calculating 

200  points  for  plotting;  an  alternate  method  is  desirable. 
Graphs  of  the  numerical  results  are  withheld  until  the  approxi- 
mate result  is  derived. 

From  Appendix  A,  the  Q function  can  be  expanded  in 
an  infinite  series  of  modified  Bessel  functions  and  using 
a and  b as  defined  in  Eqs  (48)  and  (49),  the  probability 
of  error  is  rewritten: 


Pe(A0  = 1/2  exp  ( ~R(^+B))  . IQ(c)  + 2£  (d)ni  (c)  (55) 

n = 1 


The  form  of  Eq  (55)  cannot  be  integrated  directly  using 
f (A0 ) however,  the  modified  Bessel  function  can  also  be 
expanded  in  a power  series  and  using  the  first  two  terms  in 
the  power  series  expansion  from  Appendix  A,  the  conditional 
probability  of  error  is  rewritten: 


?e|A0  = 1/2  exp  (~R(2*g))  • |2exP  (f>  * dc(exp(|)-l) 


-1  - 


(56) 


The  expression  in  Eq  (56)  can  be  averaged  term-by-term 
using  the  circular  moment  properties  of  f(A0)  (Lindsey,  1973: 
36) , resulting  in: 


Pe  = 1/2  exp  (Z^i |±£I) 


2 exp  (Ei^iiL)  *0  (a-^L.) 


0 (a) 


+ I 


(1+B  ) exp  (R(*+g)) 


+ 2 /B  -exp  (R<i±!> 


T^rr 


- d+e) 


- 2/B  Xi  (a)  - 1 - R2 


(d^)2  - 2B(1+M^_j] 

' 10(a)' 


(57) 


The  resultant  expression  is  rather  cumbersome;  however, 
the  probability  of  error  can  easily  be  plotted  using  a com- 
puter plotting  routine  and  Eq  (57).  The  range  that  Eq  (57) 
gives  valid  results  is  determined  by  comparing  the  numerical 
results  of  Eq  (54)  to  the  approximate  results  of  Eq  (57)  . The 
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results  for  various  values  of  variance  of  A0  with  6 =1 

are  given  in  Figure  8 thru  12.  From  the  results  of  Figures 
11  and  12,  the  approximation  is  within  .5  dB  of  the  numeri- 

_ 7 

cal  results  for  an  R necessary  for  P_  > 10  . The  approx- 

hi 

imation  is  of  little  use  when  a 2 - .1  as  it  is  accurate 

A0 

_ 2 

only  when  P„  >10  as  can  be  seen  in  Figure  8.  When  oA  2 = 
hi  A 0 

.05,  the  approximation  is  useful  from  R's  necessary  for 

_ 3 

P£  >10  . In  Figure  10  the  approximation  is  useful  for  R's 

_ 4 

necessary  for  P£  >10  . Table  I summarizes  these  results 

and  the  approximate  values  of  R for  a given  P with  the 
* indicating  those  values  that  the  approximation  is  within 
10%  of  the  numerical  results.  The  range  for  which  Eq  (57) 
is  useful  within  10%  of  numerical  results  is  0<  a 2 < .033 

— A 0 — 

and  R<10dB. 

Table  I . Required  R in  dB  for  a Given  Variance  and  Prob- 
ability of  Error  ( g=T) 

Numerical  Results 


PE 

°2A0  " -1 

.05 

.033 

.02 

.01 

_ 3 

10 

9.75 

8.54 

8.29 

8.11 

7.93 

-5 

10 

15 

11.78 

11.05 

10.70 

10.48 

PE 

Approximate 

Results 

10"3 

8.80 

8.36* 

8.22* 

8.10* 

8.02* 

in 

1 

O 

H 

11.35 

10.88* 

10.72* 

10.58* 

10.46* 

Benchmark  Receivers 

PE 

Coherent 

(e=“'  aA0=O) 

DPSK  ( 

2 

a =0) 

A0 

10“3 

6.79 

7.93 

10"  5 

9.58 

10.34 
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♦These  values  of  R are  within  10%  of  the  numerical  values. 

From  Table  I,  the  approximation  for  Pe  is  very  good 
when  cr,  2 < .05  radians2  (a4Q  < 12.8°)  and  R <9dB.  The 
effects  of  increasing  8 cannot  be  demonstrated  using  Eq  (57) 
as  the  approximation  is  useful  in  an  ever-decreasing  value  of 
R.  The  numerical  techniques,  however,  can  be  used  to  deter- 
mine the  effect  of  increasing  8 . Figures  13  and  14  show  the 
numerical  results  for  two  values  of  aA02  and  a few  values 
of  8 . In  the  limit,  however,  as  8-*-°°  the  conditional 
probability  of  error  reduces  to  (Stein,  1964:44): 

P , = 1/2  erfc  (/R  cos  A0)  (58) 

6 | A 6 

Averaging  Eq  (58)  results  in  an  infinite  series  expansion 
that  is  evaluated  for  various  values  of  a in  other  texts 
(Lindsey  and  Simon,  1973:313-316)  using  f(A0)  as  given  by 
Eq  (51)  . Figure  15  shows  the  results  of  integrating  Eq  (58) 
over  f ( A0 ) for  the  same  values  of  a as  given  in  Figures  8 
and  9.  A comparison  of  Figure  13  with  Figure  15  for  aA92  = •! 
and  8 = 10  indicates  that  there  is  little  additional  improve- 
ment in  performance  when  8 is  increased  from  10  to  infinity. 
This  result  shows  that  when  aAQ2>.l  the  variance  is  the  more 
sensitive  performance  parameter.  The  averaging  of  Eq  (58) 
results  in  an  irreducible  error  performance;  for  example, 
when  o is  10  or  less  (cja02  *.1  rad2)  the  receiver  will  not 
perform  better  than  10  5 errors  per  bits  for  R as  high  as 
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Figure  8.  Numerical  vs  Approximation  Probability  of 


Error 


A0 


.1 
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Figure  12.  Numerical  vs  Approximate  Probability  of  Error 


Probability  of  bit  error  (log 


Energy  per-bit-to  noise  spectral  density  (R  in  dB) 


¥ 

i 

Figure  Error  performance  with  o*e  = .02  and 

various  values  of  B 
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I 

> 1 

20dB  or  greater.  The  impact  of  these  results  will  be  dis- 
cussed in  the  final  chapter  after  some  special  cases  of 
the  generic  receiver  are  examined. 

Special  Configuration  of  the  Generic  Receiver.  As 
previously  mentioned,  if  the  filters  in  the  generic  receiver 
are  bypassed  and  the  data  feedback  eliminated,  the  resultant 
receiver  is  the  DPSK  receiver  of  Figure  16,  where  Figure  16 
is  a redrawing  of  Figure  5 relabeled  in  terms  of  the  compo- 
nents of  S and  U . The  vectors  S and  U for  the 
o o 

DPSK  receiver  are  redefined  as 


U(k)  = 


S (k)  = 
o 


(59) 


(60) 


where 

I(k-l)  = i(k)  of  the  generic  receiver  and  Q(k-l)  = 'S(k) 
of  the  generic  receiver  for  this  special  DPSK  case.  The  (k-1) 
indicates  that  the  k*"*1  data  bit  was  delayed  T^  seconds. 


The  optimum  decision  from  Eq  (43)  is: 

l 

> 

I (k)  I (k-1)  + Q (k)  Q (k-1) o 0 (61) 

The  performance  given  A0  is  given  by  Eq  (44)  or  (45) 
where  0=1  in  a and  b so  that: 
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a = R ( 1-cosA  0 ) 


(62) 


b = R (1+COSA0 ) 


(63) 


The  results  for  the  DPSK  receiver  when  phase  varies  dur- 
ing 2 seconds  is  given  by  Figures  8 thru  12  for  various 

values  of  °A02  • From  Appendix  A when  a = °°  or  a A02  = 0 
Eq  (44)  and  (45)  reduce  to  the  DPSK  benchmark  results  of  Eq 
(40),  as  would  be  expected. 

A possible  modification  of  the  DPSK  receiver  would  be  to 
estimate  the  value  of  the  received  vector  SQ  when  a binary 
one  is  sent  or  180°  if  binary  zero  is  sent. 

The  optimal  estimator  is  the  maximum  likelihood  (ML) 
estimator.  The  ML  processor  is  given  in  Figure  17.  The  deriva- 
tion for  the  ML  estimate  of  cos0 (k)  and  sin0 (k)  is  rather 
detailed,  and  is  derived  in  Appendix  C.  The  assumptions  made 
in  deriving  the  ML  estimator  are  that  the  0 (k)  phase  measure- 
ments are  independent  and  there  are  no  a priori  statistics  on 


0 (k)  . 

The  form  of  SQ  remains  unchanged,  however,  U is  given 


by: 


U = 


I cos© (k-1) ‘ 
sinMlo-1)  . 


(64) 


where 
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A 

X is  the  ML  estimate  of  X 

The  ML  estimator  of  Appendix  C merely  normalizes  the 
U vector  to  a unit  vector  and  the  detection  scheme  is  equiva- 
lent to  the  DPSK  receiver.  The  error  performance  is  also 
equivalent  to  the  DPSK  receiver.  The  ML  estimator  proceeded 
by  the  data  feedback  provides  a normalized  output  suitable 
for  further  processing  such  as  predicting  the  k*"*1  A0  from 
the  preceeding  values.  As  previously  mentioned,  the  ML  esti- 
mator assumed  there  were  no  apriori  statistics  on  A0  and 
the  0 (k)  measurements  were  uncorrelated.  The  phase  process 
0 (t)  is  however  well  defined  and  the  correlation  given  by 
Eq  (9)  it  is  then  conceivable  that  useful  information  was 
ignored  in  using  the  ML  estimator.  Perhaps  a more  suitable 
estimator  is  a minimum  mean  squared  error  (MMSE)  or  a maximum 
posteriori  (MAP)  or  a predictor  such  as  a Kalman  filter.  Such 
predictors  and  estimators  are  not  the  subject  of  this  paper, 
however,  they  could  possibly  reduce  a2AQ  which  from  Figures 
8 thru  12  clearly  improves  performance  of  the  receiver. 
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Chapter  IV.  Conclusions  and  Recommendations 


The  significant  results  and  recommendations  for  fur- 
ther study  are  addressed  in  this  chapter. 

Conclusion 

The  advantage  of  examining  the  generic  receiver  is 
evident  in  the  limiting  cases  and  by  varying  aA02  and 
3 . The  graphs  of  performance.  Figures  8 thru  12,  show 

that  °Ag2  is  the  most  sensitive  parameter.  When  8 is 
varied  from  1 to  » , the  performance  improves  slightly  in 
the  threshold  region,  whereas  when  a.  2 is  reduced  from 
.1  rad2  to  .05  rad2  , the  performance  improved  rapidly. 
The  results  indicate  that  the  design  effort  in  the  generic 
receiver  should  be  concentrated  on  predicting/estimating 
the  phase  quadrature  to  reduce  aA62  • The  apriori  model 
for  the  differential  phase  A0  is  well  developed  there- 
fore, it  is  conceivable  that  a Kalman  filter,  MAP,  or 

MMSE  could  reduce  a.  2 . 

A o 

The  requirement  that  CCDs  operate  at  baseband  indi- 
cates that  the  generic  baseband  receiver  is  a useful 
approach  to  pursue.  The  CCDs  are  an  excellent  choice  for 
code  correlators  because  of  their  flexibility,  however, 
they  operate  in  the  baseband  region  and  require  processing 
of  the  incoming  signal  at  baseband.  Because  of  this  base- 
band processing  restriction,  the  generic  baseband  receiver 
is  not  the  optimal  processor.  Additionally,  the  perform- 
ance is  completely  defined  by  oA02  and  6 , thus  simpli- 

fying performance  analysis. 
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Recommendations 

The  following  topics  for  future  study  are  not  com- 
plete, however,  they  represent  the  logical  extension  of 
this  thesis. 

The  statistical  model  for  the  differential  phase 
process  is  well  defined  in  Chapter  II.  The  statistics  of 
the  phase  process  could  be  exploited  in  the  design  of  a 
Kalman  filter,  MAP,  or  MMSE  predictor/estimator  for  the 
quardratures  of  this  receiver.  The  effort  in  this  approach 
would  be  to  reduce  a.  2 and  increase  B through  predic- 
tion  or  estimation  on  the  quadratures.  An  extension  to 
this  investigation  is  relaxing  the  restrictions  on  Tc 
over  Tfa  and  possibly  introducing  sub-bit  processing. 

The  next  logical  step  in  investigation  is  the  intro- 
duction of  a more  sophisticated  doppler  model.  The  receiver 
could  possibly  be  aided  by  inertial-navigation  system  (INS) 
or  other  information  external  to  the  receiver  to  improve 
performance  with  the  more  sophisticated  phase  model. 

From  the  recommendations,  it  is  obvious  that  the 
analysis  is  not  complete.  This  paper  set  up  the  sampled 
data  model  and  developed  the  baseband  processor.  The  per- 
formance did  not  result  in  a simple  expression.  As  was 
mentioned  in  the  introduction,  the  price  for  increased 
resistance  to  jamming  is  an  increase  in  receiver  complexity. 
Unfortunately,  this  usually  translates  into  increased  diffi- 
culty in  pinpointing  error  performance. 
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Appendix  A 


The  Q Function  and  Modified  Bessel  Function 
This  appendix  contains  some  useful  relations  for  the 
Q function  and  modified  Bessel  function  (Shwartz,  1966: 
585-589;  and  Van  Trees,  1968:394-396).  The  definition  of 
the  Q function  is  given  by: 


Q(a,b)  = 


■/ 


exp 


I.  (ax)x  dx 


(A-l) 


where 


IQ  is  the  modified  Bessel  function  of  order  zero. 

For  reference,  the  modified  Bessel  function  of  order  n 
is  given  by: 

I (x)  = /”2ir  exp  (+jn6 ) exp  (xcose ) de  (A-2) 

n Z1T  Jo  ~ 

The  limiting  properties  of  the  Q function  are: 


Q (o,b)  = exp 

(2_) 

(A-3) 

Q (a,o)  = 1 

(A-4) 

Some  useful  approximations  of  the  Q function  are: 

Q(a,b)  = 1/2  erfc  b>>l,  b>>b-a  (A-5) 

/2 
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The  series  expansions  for  the  Q function  are  given  by: 


Q(a,b)  = exp  (~a  *b  ) E 

n = o 

(g)n 

Vab> 

a<b 

(A—  6 ) 

Q(a,b)  = i-exp(-^i^i)  £ 

n = 1 

(|)n 

Vab) 

b>a 

(A-7) 

Three  symetric  and  antisymetric  relations  for  the  Q 
function  are  given  by: 


1 + Q(a,b)  - Q(b,a)  = f “ exp(-y)In  (2ab£_)dy 

b2+a2  J2  2 0 a2+b2 

a + b 



b>a>o  (A-8) 


1 + Q (a,b)  - Q(b,a) 


erfc  (^=^) 
✓2 


b>>l,a>>l,b>>b-a>0 


(A— 9 ) 


(A-10) 


The  series  expansion  for  the  modified  Bessel  function 
is  given  by: 


In<x> 


00 

k=o 


x n+2k 
1 <*) 
kf  (n+k)  1 2 


(A-ll) 
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{ 


I (X)  * 


n 


/2ttx 


1 - 


4n2-l 

8x 


x>>  1 


(A-13) 


Appendix  B 


( 


A Recursive  Method  of  Computing  the  Q function 

The  Q function  can  be  computed  recursively  to  a given 
accuracy  by  using  the  following  equations  (Brennan,  1965: 
312-313) , 


Q(a,b)  = 1 -21  g (n)  *k (n) 
n=o 


where : 


and: 


k<„)  = |i  klszll 


with: 


g(o)  ■=  e~b2/2 


and: 


k(o)  - e 


-^2 


az/2 


(B-l) 


g(n)  - g(n-l)  - ^ (|i)n  e-b2/2  (B-2) 


(B-3) 


(B-4 ) 


(B-5) 


After  N iterations,  the  error  in  the  remaining  terms 
is  given  by: 
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i- 


1 


Rv,  = g(n)  *k(n) 
N 


(B-6 ) 


A useful  upper  bound  on  Rjj  is  given  by: 


Rjj  < (k(N)  • g ( N ) / ( 1 — ( — — — ) ) 
" 4N2 


( B— 7 ) 


for: 


N > 


ab 


(B-8) 


The  summation  in  Eq  ( B— 1 ) is  initialized  using 
Eqs  (B-4)  and  (B-5) , then  Q is  computed  recursively  until 
the  remainder  criteria  in  Eq  (B-7)  is  met  or  exceeded.  The 
series  converges  rapidly  when  n>N  (Brennan,  1965:313). 
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Derivation  of  Maximum  Likelihood  Estimate  of  Cos  9 and  Sin  9 

The  receivers  in  Chapter  III  all  perform  some  form  of 
the  dot  product  to  recombine  the  quadrature  signals  from 
the  output  of  the  correlators.  One  method  of  obtaining  the 
projection  of  the  received  signal  vector  So  onto  U is 
to  form  the  scalar  product.  The  scalar  product  of  two  vec- 
tors So  and  U is: 


So 


SQ  cos  ^ 


where 


|U|>0  (C-l) 


x is  the  vector  x 
| x | is  the  magnitude  of  the  vector 
x • y is  the  dot  product  of  the  vector  x and  y 
0 is  the  smallest  angle  between  So  and  U 

The  quantity  — — is  a unit  vector  in  the  direction 

|U  | ^ 

of  U . The  components  of  b are  defined  as  in  Figure 
C-l  and  the  components  of  the  unit  vector  are: 

= coseT  + sine  J (C-2) 

|U  | 

where 


0 is  the  angle  between  U and  Sq  direction 
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Figure  C-l . Dot  Product  of  Vectors  S and  U 

o 


So  + — = S cosel  + Sn  sinej  (C-3) 

|U|  0!  02 


where 


So  1 

is 

the 

component 

of 

So 

in 

the 

T 

direction 

Sq  2 

is 

the 

component 

of 

So 

in 

the 

3 

direction 

The  components  of  So  can  be  rewritten  as: 

S = | SqI  cos'?  (C-4) 


57 


S02  = I S0  | sin1?  (C-5) 

i 

Using  Eqs  (C-4)  and  (C-5)  in  Eq  (C-3) : 

( C— 6 ) 

which  reduces  to: 

S0  . — = IS  I cos(V-e)  (C-7) 

|U|  I 

But,  y-0  = 0 from  Figure  C-l  and  Eq  (C-7)  is  rewritten: 

U 

s0  • |u|  = lsol  cos  * (c_8) 

The  results  in  Eq  (C-8)  are  the  motivation  for  the 

determination  of  the  cose  and  sine  components  of  the  vector 

U . In  Figure  C-l,  the  vector  So  represents  the  received 

signal  r(t)  decomposed  into  an  orthogonal  vector  rc  and 

r dotted  with  the  signal  vector  sg  • The  U vector  is 
s 

defined  in  a similar  manner  except  it  dots  the  r and  r 

c s 

vector  with  the  u vector.  The  closer  the  estimate  of  cose 
and  sine  is  to  cose  and  sinf  the  more  closely  the 
receiver  approaches  coherent  operation  with  e = f the 
coherent  result. 

The  maximum  likelihood  (ML)  estimate  of  cose  and  sine 
given  the  statistics  of  U is  the  best  estimate  of  cose 
and  sine  (Van  Trees,  1968:160).  The  ML  estimate  is  found 

( 

by  solving  the  equation  for  the  sufficient  condition: 
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dlnf (x ly) 
dy 


0 


( C— 9 ) 


*"*ML 

where 


d is  the  partial  derivative  operator 
ln(x)  is  the  natural  log  of  x 


The  output  of  the  correlators  is  I(k)  and  Q(k) 
as  given  in  Chapter  II.  Using  the  reference  in  Figure 
C-l,  I(k)  is  in  the  I direction  and  Q(k)  is  in  the 
7 direction.  The  conditional  density  of  I (k)  and  Q(k) 
given  8 and  that  the  noise  is  independent  between 
samples  is: 


N 


f ( i r Q I e ) = n exp 

k=l  z a 


N 


-1/2  ^ -ANTccos9 (k>I  2 * * * * * * 


’k=l 


-1/2  (Q(k)  -ANTcSi*10  00  ) 2 


(C-10) 


where 


I is  an  NX1  vector 

Q is  an  NX1  vector 

<f  is  an  NX1  vector 

N 

n is  an  N-fold  multiplication  of  each  term  from 

k=l 

k = 1 to  N 

a2  is  the  variance  of  the  Gaussian  variables 
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The  data  term  d(k)  in  I(k)  is  assumed  to  be 
di (k)  so  that  Eq  (C-10)  is  a valid  density  to  use  in 
Eq  (C-9) . If  the  phase  is  assumed  to  be  independent 
between  samples,  the  ML  estimate  for  cose (k)  is  found 
by  solving: 


dlnf (I,Q|e ) _ 

dcose (k) 

cose  (k)  = cose(k)„_ 

NJ-i 


(C-ll) 


The  results  in  Eqs  (C-14)  and  (C-15)  seem  intuitive 
from  the  operation  in  Eq  (C-2) . Additionally,  the  estimates 
are  the  ML  estimates  for  cos6 (k)  and  sine (k)  . The  assump- 

tions for  Eqs  (C-14)  and  (C-15)  to  be  valid  are  that  the  phase 
0 (k)  is  constant  over  a bit  interval  T^  . Also,  the  phase 
0 (k)  is  independent  of  the  previous  phase  sample. 

The  Cramer-Rao  (C-R)  bound  of  the  estimates  of  the 
cose  and  sine  give  an  indication  of  the  variance  of  the 
lower  bound  on  the  estimate.  The  C-R  bound  is  given  by  (Van  Trees, 
1968:437): 


d2f  (x(y)1 


(C-15) 


where 


E |xj  is  the  expected  value  of  x 
d2  is  the  second  partial  derivative 


The  first  partial  derivative  of  cos0  is  the  right  side 
of  Eq  (C-13)  so  taking  the  partial  derivative  of  Eq  (C-13)  and 
letting  0 = cos0 (k)  and  sin0  (k)  = /^_g 2 yields 


o<2f  (I, Ole) 

d cos0 (k) 2 


Q (k) ANT 
” c 

a2  (1-0)2 


1-02 


1-02 


(C-16) 


Taking  the  expected  value  of  Eq  (C-17)  and  inverting: 


cos0 (k) , 


Nq  sin20 (k) 


(C-17) 
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where 


2^  is  given  as  the  energy  per  bit  from  Chapter  II, 


The  C-R  bounds  can  be  rewritten: 


C-R^  „7b\  - R sin2  0 (k) 

cose (k) ml 


(C-19) 


C_Rslne7k)‘ML  = R cos2  6 (k) 


(C-20) 


The  C-R  bounds  for  cos0(k)„T  and  sin0(k),.T 

ML  ML 


function  of  0 (k) 


are  a 

Unfortunately,  as  in  the  case  of  0 (k)  = 0 


radians,  the  C-R  bound  for  cos0(k),,T  is  zero,  which  underbounds 

nil 

all  estimates.  The  C-R  bound  in  this  case  does  not  give  a tight 
bound  on  the  variance  of  the  expected  performance.  This  fact 
does  not  exclude  the  use  of  Eqs  (C-14)  and  (C-15)  as  the  optimal 
estimates  it  is  just  that  tight  bounds  on  performance  cannot  be 
determined  using  the  C-R  bounding  technique. 
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for  specific  cases  and  is  useful  for  small  values  of  variance  and  low  signal-to- 
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noise  ratios.  In  the  limit  when  the  variance  approaches  zero,  the  receiver 
performance  expression  is  equivalent  to  the  differentially  coherent  phase 
shift  keyed  (DPSK)  receiver.  If  the  reference  is  perfect  in  addition  to  the 
zero  phase  variance  the  performance  is  the  same  as  a phase  coherent  receiver 
with  a perfect  reference. 
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